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Tank-treading, tumbling and trembling are different types of the vesicle behavior in an external 
flow. We derive a dynamical equation for nearly spherical vesicles enabling to establish a phase 
diagram of the system predicting the regimes. The diagram is drawn in terms of two dimensionless 
parameters depending on the vesicle excess area, fluid viscosities, membrane viscosity and bending 
modulus, strength of the flow, and ratio of the elongational and rotational components of the flow. 
The tank-treading to tumbling transition occurs via a saddle-node bifurcation whereas the tank- 
treading to trembling transition occurs via a Hopf bifurcation. We establish a critical slowing near 
the merging point of the transition lines. 
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I. INTRODUCTION 

Vesicles, which are closed membranes separating two 
regions filled up by generally different fluids, have very 
much in common with biological cells and have a num- 
ber of applications in pharmaceutics. The truly non- 
equilibrium problem of describing the vesicle dynamics 
in external flows attracted lots of attention both from 
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1231 . I24I . I25j . In experiments different regimes of vesicle 
dynamics were observed, usually in share flows. In the 
tank-treading regime a vesicle shape is stationary. In the 
tumbling regime a vesicle experiences periodic flipping 
in the shear plane. Trembling, experimentally discovered 
in the work [7|, is an intermediate regime between tank- 
treading and tumbling, in the regime a vesicle trembles 
around the flow direction. Theoretical studies [l^l and 
[20l | have predicted similar vesicle behavior and called it 
vacillating-breathing and swinging respectively. The dif- 
ferent types of the motion are illustrated in Fig. [1] 

In this paper we focus on constructing the phase di- 
agram, which would predict what kind of regime corre- 
sponds to a given set of external parameters, such as 
vesicle non-sphericity, viscosity contrast, strength of ex- 
ternal flow. Although this problem was addressed pre- 
viously by different authors, it is still far from being 
closed. As long as no analytic solution of this problem 
exists, theoretical studies were based either on numerical 
simulations or on some approximations allowing analyt- 
ical treatment. Numerical investigations of this problem 
involved several different computational schemes, includ- 
ing boundary element method [alSl , mesoscopic particle- 
based approximation l^j. 12Jl3 , [ij] , and an advected 
field approach [H, These approaches have 

shown qualitative agreement with experiment. However, 
they did not solve the problem of constructing the vesi- 
cle phase diagram completely. Analytical studies of the 
problem can be divided in two major classes. In the first 
one, see Refs. [ll|, El, HSl, phenomenological models of 
vesicle dynamics based on the classical work of Keller 
and Skallak [2l[ were proposed and proved themselves to 




FIG. 1: Vesicle projections to the shear plane in the tank- 
treading, trembling and tumbling regimes. 



be rather efficient in explaining particular experiments. 
However, they cannot cover wide range of parameters and 
cannot, consequently, be used for constructing a phase di- 
agram for the vesicle dynamics. In the second series of 
works, Refs. [22i, i23„ 2A, 2Mi the studies are focused on 
nearly spherical (quasi-spherical) vesicles whose shape is 
close to spherical one and can be parameterized by an 
expansion over spherical harmonics. 

We propose a natural extension of the theory devel- 
oped for quasi-spherical vesicles by accounting higher- 
order expansion terms. A perturbation scheme around 
the Lamb solution for spherical body in an external flow 
[3^ allows one to derive the dynamic equations for the 
shape and the orientation of a vesicle and to investigate 
them analytically. We show that the high-order in the 
vesicle distortion terms produce a qualitative change in 
the phase diagram and make it significantly more com- 
plicated. The resulting diagram contains all three types 
of vesicle behavior which were observed experimentally. 
In our work, we analyze how the vesicle dynamics de- 
pends on different control parameters, such as viscosity 
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contrast, vesicle excess area, internal membrane viscos- 
ity, strength of the flow, and ratio of the elongational and 
rotational components of the flow. We analyze also the 
vesicle orientation in the tank-treading regime. Some of 
the results, derived in this paper were already reported in 
p8l |. Here we present significantly more detailed deriva- 
tions and analyze several aspects of the problem which 
were not discussed in (isj . 

Speaking about membranes we have in mind lipid bi- 
laycrs, the simplest type of the biological membranes. 
Physical properties of such objects have been extensively 
studied, both experimentally and theoretically (see e.g. 
the books [H,!!!!!! and the reviews [H, [13, [Sfl ) • There 
are several features of the membranes which are impor- 
tant for our analysis. First, we assume that the mem- 
brane is in a fluid (isotropic) state (is a 2d liquid), which 
is typical of lipid bilaycrs under normal conditions. Sec- 
ond, we assume that the vesicle has an excess area, that 
enables one to treat the membrane as incompressible. 
Third, we assume that the membrane is impermeable to 
the surrounding fluids, the condition is usually well sat- 
isfied in experiment. Finally, we take into account the 
membrane internal viscosity, which can play an essential 
role, say, in the vicinity of the lipid-bilayer melting point 

The structure of our paper is as follows. In Section 
|TT]wc expose basic theoretical facts concerning physics of 
membranes. A special attention is paid to dynamic prop- 
erties of the membranes which can be treated as moving 
interfaces immersed in a 3d fluid. In Section IlIII we for- 
mulate peculiarities of nearly spherical vesicles, analyze 
their equilibrium properties, and derive phenomenologi- 
cal equation which describe their dynamics in weak exter- 
nal flows. In Section Hv] we establish a dynamic equation 
for the vesicle evolution obtained in the framework of ex- 
pansion over deviations from the ideal spherical shape. 
We introduce dimensionlcss parameters controlling the 
vesicle behavior. In Section |V] we consider the case of 
planar external velocity field. In this case the dynamic 
equations can be essentially simplifled enabling one to 
establish the phase diagram of the system containing 
domains corresponding to tank-treading, trembling and 
tumbling. In Section IVII we analyze limit cases of weak 
and strong external flows where it is possible to study 
the vesicle behavior in details. In Conclusion we discuss 
some outcomes of our work and its possible extensions. 
Some technical details of the perturbation expansion are 
presented in Appendix. 



II. BASIC RELATIONS 

We are interested in processes which take place at 
scales of the order of the vesicle size which are assumed 
to be much larger than the membrane thickness. This 
assumption is well justified for giant vesicles, usually ex- 
amined in experiment. In the main approximation, the 
membrane can be treated as infinitcsimally thin, that is 



as a 2d object (film) immersed into a 3d fluid. In this 
case the vesicle is characterized by its geometrical shape. 
In other words, in this limit the vesicles can be consid- 
ered as active interfaces separating different pieces of the 
fluid. 

The two membrane properties, its incompressibility 
and impermeability, imply that both the vesicle volume V 
and its surface area A are conserved, provided the vesicle 
has an excess area. The latter can be characterized by a 
dimensionlcss factor A, which is traditionally introduced 
as 

A^{ATT + A)rl V^iTT4/3, (1) 

where ro is a vesicle "radius" determined by its volume. 
The excess area factor is non-negative, A > 0, and the 
minimal value A = corresponds to the ideal spherical 
geometry. The nearly spherical (quasi-spherical) vesicles 
are characterized by the condition A <C 1. 

The energy of the membrane is determined by its bend- 
ing distortions and can be written as the following surface 
integral [H [H, [H, 113 

J-(b) = J dA {k.H^/2 + RK) , (2) 

over the membrane position. Here n and R are bending 
modules, H and K are the mean and Gaussian curva- 
tures, respectively. They are related to the local curva- 
ture radii of the membrane i?i and R2 as 

= K^R^^R^^. (3) 

In accordance with the Gauss-Bonnet theorem, the sec- 
ond term in the right-hand side of Eq. ^ (with the 
factor R) is invariant under smooth deformations of the 
membrane shape. Therefore the term is irrelevant pro- 
vided the vesicle topology is fixed. 

Besides the bending energy ^ the membrane is char- 
acterized by its surface tension a. One should stress that 
for vesicles with excess area the surface tension is an aux- 
iliary variable adjusting to other quantities characteriz- 
ing the membrane and the surrounding flow to ensure the 
membrane incompressibility. 

A. Flow near the vesicle 

In this paper we consider the situation where both flu- 
ids, outside and inside the vesicle, are Newtonian. Fur- 
thermore, we assume that the Reynolds number associ- 
ated with the vesicle dynamics is vanishingly smalL which 
is the case in real experiments [l|, 0, H, m H, B ■ Un- 
der these assumptions the fluids can be described by the 
Stokes equation 

gdtv = T]S/^v-\/P , (4) 

where P is pressure, v is the fluid velocity, g is its mass 
density, and rj is its dynamic viscosity. The equation 
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(m has to be supplemented by the incompressibihty con- 
dition Wv = 0, which leads to the Laplace equation 
V^P = for the pressure. 

We divide the flow near the vesicle into two parts: an 
external flow, which would be observed in the fluid in 
the absence of the vesicle, and an induced flow, which 
is excited as a result of the vesicle reaction to the exter- 
nal flow. The external flow is assumed to be stationary, 
its velocity is designated as V. One should remember 
that the vesicle is advected by the flow, and therefore 
the above assumption should be valid in the Lagrangian 
reference frame attached to the vesicle. Below, we ne- 
glect the term with the time derivative in Eq. (|4|) since 
the characteristic time scale associated with the vesicle 
dynamics is assumed to be large compared to the viscous 
relaxation time Qr^/rj. 

We assume that the characteristic spatial scale of the 
external flow is much larger than the vesicle size. In 
this case the external velocity V near the vesicle can be 
approximated by a linear profile, determined by a deriva- 
tive matrix d^Vi. The incompressibihty condition implies 
that the matrix dkVi is traceless. Generally, the external 
flow has two contributions, elongational and rotational: 

dkVi ^ Sik - eikjUJj , (5) 

where s is the strain matrix (symmetric part of the ma- 
trix dkVi) and a; is the angular velocity vector. The 
strain can be characterized by its strength s, defined as 

= (l/2)tr s^. Note that for a shear flow s = = 7/2, 
where 7 is the shear rate. 

The fluids inside and outside the vesicle are assumed to 
be different. We use the designations 77 for the external 
fluid viscosity, the viscosity of the internal fluid is desig- 
nated as 77. An important parameter which controls the 
tank-treading to tumbling transition is the viscosity con- 
trast 77/77. The limit where the viscosity contrast tends 
to infinity corresponds to a solid body behavior of the 
vesicle, that was concerned by Jeffery |46j . 

The membrane moves together with the fluid that is 
the velocity field v is continuous on the membrane and 
the field v determines the membrane velocity as well as 
the fluid velocity. For relatively slow processes, we are 
investigating, the membrane can be treated as locally 
incompressible, that leads to the condition 

d:^Vi = 0, where d:^ — S^dk , (6) 

to be satisfied on the membrane. Here 5^ is the projector 
to the membrane, it can be written as 5^ = Sik — Idk 
where I is the unit vector normal to the membrane. The 
3d incompressibihty condition = together with Eq. 
^ leads to the relation UlkdiVk = 0, to be satisfied at 
both sides of the membrane. 

B. Membrane stress 

The membrane reaction is characterized by its surface 
stress tensor T^^ . There are three contributions to the 



tensor related to the bending energy ([2]). to the surface 
tension of the membrane, and to the internal membrane 
viscosity: 

= ' - ^4 - QSi^St + d^v,) , (7) 

where a is the surface tension coefficient and Q is the 
membrane (2d) dynamic viscosity. Note that the surface 
tension a plays an auxiliary role being adjusted to other 
stresses to ensure the local membrane incompressibihty. 
An expression for the bending contribution to the surface 
stress tensor was found in the work [i^ (see also the book 
[3^). It can be written as 

t'i^ = « (- ^ff '4 + Hdih kdiH^ , (8) 

where = VZ is the membrane mean curvature and 
is defined by Eq. 

The surface force / (force per unit area) associated 

fs) 

with the membrane stress tensor can be calculated 

as fi = —d^T^^\ Then one obtains from the expressions 
(|7l8p three contributions to the surface force 

/ - f^^^ + + , (9) 

where 

/f = k[H - 2K) + A^H] I, , (10) 

fl"^ ^-Hok + dtcT, (11) 

/f ) = C - HLdtvn - 2hid^l,)d^vn] . (12) 

Here, again, H and K are the mean curvature and the 
Gaussian curvature of the membrane, and A-*- is the 
Laplace-Beltrami operator. A-'- = d^d^, associated with 
the mebrane. Note that the expression for the force (fTO|) 
can also be derived by calculating the variation of the 
bending energy Q due to infinitesimal membrane defor- 
mations [sst . 

The surface force / is compensated by the momen- 
tum flux from the surrounding medium to the membrane. 
This flux consists of two parts, related to the fluid pres- 
sure and the fluid viscosity. As a result of the balance, 
we find the following relations 

- k[H{H^/2 - 2K) + A-^F] + aH 

+2Cdilndiv, ^ P„, ~ Pouu (13) 

= h [vidtVk + dkVi)in - vi^iVk + dkVi)out], (14) 

for the normal and tangential to the membrane compo- 
nents of the force. Here, we assumed that the unit vector 
I is directed outwards the vesicle and the subscripts in 
and out label regions inside and outside the vesicle, re- 
spectively. Thus, Pin — Pout is the pressure difference 
between the inner and outer regions that is the pres- 
sure jump on the membrane. Note that a fluid viscous 
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contribution is absent in Eq. (|13p due to the condition 
liljdiVj = 0, following from the membrane incompress- 
ibility (see above). 

To find the velocity field at a given membrane shape 
one should solve the stationary Stokes equation 7]\/^v — 
VP (inside and outside the vesicle) with the boundary 
conditions (|6ll3ll4p on the membrane. An additional 
boundary condition reads that v —^ V far away from the 
membrane. Note that due to linearity of the equations 
and the boundary conditions for the velocity a solution 
of the equations can be written as a sum 

V = v^'^ + v'^'^l (15) 

where t?*-*^ is proportional to the gradient of the external 
flow (O, and v'-'^') is proportional to the bending modulus 
K. Of course, t;'-*' and v'-'^'> are complicated functions of 
the vesicle shape. 



C. Membrane shape parametrization 



displacement u{9,(p). Second, one should use the kine- 
matic relation ([TO]) . This procedure results in a closed 
non-linear equation for u. Note that due to the prop- 
erty (jlSp the expression is a sum of two terms, which are 
proportional to the external flow gradient ([5]) and to the 
bending modulus n. Of course, both terms are non-linear 
in u. 



III. PERTURBATION EXPANSION: WEAK 
FLOWS 

Below, we consider nearly spherical (quasispherical) 
vesicles, that is the excess area parameter A, introduced 
by Eq. ([1]), is considered to be small. In this case the 
dimensionless displacement u is small and it is possible 
to develop an expansion over u. This perturbation series 
is a basis for subsequent consideration. 

It is natural to represent the function u(6, ip) as a sum 
over spherical harmonics: 



Below, we use a particular parametrization of the vesi- 
cle shape which is 

r = ro[l + w(0,(^)], (16) 

where tq is determined by the relation ([T|). We use 
the spherical coordinate system with variables r, 9, (p and 
with the origin in the center of the vesicle. The dimen- 
sionless radial displacement u characterizes the devia- 
tions of the membrane shape from the spherical one. 

There are constraints imposed on the function m(0, tp) 
due to the conditions ([T]). In terms of the displacement 
u the conditions can be rewritten as 

j dip dO sin e{u + u^ + = , (17) 

A = J dipde sme^^{l + u) (18) 

[(1 + uf + {du/def + sin-2 9{du/dipf] - l} . 

The relations (|17ll8p are formally exact. However, they 
can be directly used only iiu{6, ip) is a single- valued func- 
tion. 

Advection of the membrane by the surrounding fluid 
implies the following kinematic relation 

dtU = —Vr - - ( VgdgU + ) . (19) 

ro r \ smO ^ ^ J 

Here Vr,V0, v^p are spherical components of the velocity v 
taken at the membrane, that is at r determined by Eq. 
(fTB)) . Again, the relation (|19p is formally exact, but can 
be directly used only if u{9, ip) is a single- valued function. 

Closed equation which describes the dynamics of the 
membrane displacement u can be derived in two steps. 
First, one should find the fluid velocity profile for a given 



^^ = ^^^i,m3^i,m(e,</')- (20) 

The homogeneous contribution to u (its zero angular har- 
monic wo,o) can be expressed via the inhomogeneous one 
(non-zero harmonics) from the relation (jl7p which re- 
flects the volume conservation. Substituting the obtained 
expression for the zero angular harmonic into Eq. psp we 
obtain an expression for A whose expansion over u starts 
from the second order term. Therefore the displacement 
u can be estimated as a/A- That justifies the expansion 
over u. 

Different angular harmonics in u play different roles. 
The zero harmonic mo,o can be excluded from the be- 
ginning, as we explained. First order harmonics Mi,,n, 
describe a shift of the vesicle without changing its shape, 
and therefore do not play any important role in the vesi- 
cle dynamics. The most essential role is played by the 
second angular harmonic which determines mainly the 
vesicle shape (at small A). As to higher harmonics, they 
relax fast in comparison with the relatively slow dynam- 
ics of the second harmonic. Therefore, the higher har- 
monics also do not play an essential role in the vesicle dy- 
namics. To avoid a misunderstanding, let us stress that 
the last assertion is valid for stationary external flows. As 
it was discovered experimentally [iot and explained theo- 
retically [soj . at some conditions (abrupt inversion of the 
external purely elongational flow) high angular harmon- 
ics are generated, the phenomenon is called wrinkling. 

The vesicle shape depends on the strength of the exter- 
nal flow. In weak flows it is close to an equilibrium one 
whereas in strong flows it is determined by the velocity 
gradient matrix ([5]). In this section we consider the first 
case. We discuss the equilibrium vesicle shape and then 
develop phenomenology for the vesicle dynamics in weak 
flows. 
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A. Equilibrium 

111 the absence of the external flow, an equilibrium 
shape of the vesicle can be found by minimization of an 
effective free energy 



(21) 



where the first term is determined by the expression ([2|) 
and Tg A is the membrane excess area expressed in terms 
of the displacement u. The Lagrange multiplier ct, related 
to a fixed value of the membrane area, coincides with 
the equilibrium value of the surface tension. The second 
Lagrange multiplier (related to the volume V) is absent in 
Eq. (j21[) since we imply that the zero angular harmonic 
in an expansion of the displacement u is expressed via 
other ones from the relation ([T7)) . Therefore, the volume 
conservation is automatically satisfied in our scheme. 

If A is small, the principal contributions to the energy 
([2]) as well as to the excess area are of the second order 
in u. It is convenient to write the contributions in terms 
of the coefficients ui^m of the expansion (|^D|) of u{6, ip) 
over the angular harmonics: 



= f E (^ + 2)(/ + !);(/- 1)1^ 



(22) 



l>2.; 



Note, that the first angular harmonic (with Z = 1) is ab- 
sent in the expansions. The reason is that it corresponds 
to a vesicle shift as a whole, which docs not change the 
energy and the area of the vesicle. As it follows from 
Eq. ([m, the free energy is minimal if only the second 
angular harmonic is excited. In this case ctq = —(Sk/t^, 
which is an equilibrium value of the surface tension. 

Note that the second order term ([22|) is degenerate in 
TO. Therefore, in order to determine the vesicle equilib- 
rium shape, one should take into account terms of higher 
order in the expansion of the effective free energy (j2ip . 
which violate the degeneracy. In the main approxima- 
tion, it is possible to keep third order in u terms, and a 
contribution to u related to the second angular harmon- 
ics. 

For a subsequent analysis, it is convenient for us to use 
the following real basis 
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(Scos^ei-i), V2 



^/T5 



2V^ 



2V^ 



sin(26') siniy9 , V4 



sin(20) cos if , 
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sin^ 6 cos(2</j) , 



"05 = sin^ 9 sin(2ip) . 



(23) 



instead of the traditional angular functions 3^2, m- The 
functions "0^ are normalized as 



dip dO sin Oip^tpv = 



(24) 



In terms of the functions ([23)) . the contribution to u re- 
lated to the second angular harmonic, can be rewritten 
as follows 



5 

E 

u=l 



(25) 



where Ui, are some real coefficients. 

Expanding the bending energy ([2]) and the excess area 
A upto the third order in u and substituting there the 
expansion (|25p one obtains 



.^(3) = (u^u^ - E^,xu^,u,ux) + ar^ A(3) , (26) 
A^^^ = 2u^M^ - 2E^^xUf_,u^ux/3 , (27) 

where summation over repeated indices is implied and we 
designated 



dip dd sin 9 4'iJ.''Pu^\ ■ 



(28) 



Components of the object are of order unity, they 

can be found from the definition ((28l) after substituting 
the expressions (|23| . 



Minimizing the free energy ([^5]) over and deter- 
mining the Lagrangian multiplier a from the condition 
A = A(3) , one obtains 



1 



6k 



140F 



(29) 



This is the correction related to the third order term in 
the expansion of the free energy. The minimum of the 
energy corresponds to a prolate uniaxial ellipsoid. If the 
principal axis of the ellipsoid is directed along the Z-axis 
its shape is determined by the expression ui = A/2, 
that is 



4V27r 



(1 - 3cos2 



(30) 



Substituting the expression ((29)) into Eq. ((26)) wc find 
that the coefficient before u^u^ term in the effective free 
energy is estimated as K\fK. It contains an extra small 
factor -s/A in comparison with the natural estimation k. 
Thus, both, second order and third order, terms in the 
free energy ((^5)) are of the same order. That gives a 
formal justification of taking the third order term into 
account despite the smallness of u. 



B. Weak external flow, phenomenology 

Here, we analyze the case of weak external fiows that 
cannot significantly distort the vesicle equilibrium shape. 
As we established in the preceding subsection, the equi- 
librium shape of a nearly spherical vesicle is the prolate 
ellipsoid possessing the uniaxial symmetry. An orienta- 
tion of such ellipsoid in space can be characterized by a 
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FIG. 2; Reference frame related to shear flow. 

unit vector n directed along the principal axis of the el- 
lipsoid. If n is directed along the Z-axis then the vesicle 
shape is determined by the expression (|5D|) . Note that 
the vectors n and —n describe the same physical state 
since the ellipsoid is invariant under inversion. 

One can formulate a phenomenological equation for the 
dynamics of n in a weak external flow: dtfii ~ DijkdkVj, 
where dkVj is the velocity gradient matrix of the exter- 
nal flow and Dijk is some tensor related to the vesicle 
orientation. Due to the symmetry n — > — n the tensor 
Dijk contains only odd powers of n. Using the relation 
= 1, that is UiDij^ = 0, we arrive to the following 
equation 

dtUi = [{ukSij - njSik)/2 
+D {nkSij/2 + njSik/2 - ninjnk)]dkVj , (31) 

containing a single dimensionless parameter D. Deriv- 
ing this equation, we exploited the fact, that the vesicle 
dynamics should be purely rotational dtn = u) x n in a 
case of an external flow djVi = —eijktOk corresponding to 
a solid rotation. The factor D in Eq. (|3ip is dependent 
on relative significance of the viscous mechanisms and 
on the excess area parameter A, the dependence will be 
established further, see Section IVll 

For an external shear flow, it is convenient to utilize 
the following parametrization of the unit vector n 

n = (cos?9cos0, cosiJsin^, sint?) . (32) 

The components here are written in the Cartesian refer- 
ence frame attached to the flow: the X-axis is directed 
along the velocity and the Z-axis is directed opposite to 
the angular velocity vector lj (see Fig. [2]for clarification). 
Substituting the expression ([32|l into Eq. (|31]) and taking 
into account that the shear velocity gradient matrix has 
the only non-zero component dyVx ~ j one obtains 

7-19*0 =(D/2)cos(2(/))- 1/2, (33) 
j-^dti9 = -{D/A) sin(2i?) sin(20) . (34) 



Note that the dynamics of the angle (j) is separated. The 
equations (|33l34p resemble equations for a single polymer 
dynamics examined in Ref. [4J|. 

The equations (|33l34p lead to either tank-treading or 
tumbling regimes of vesicle motion. For \D\ > 1 the 
tank-treading regime is realized, with a steady tilt angle 
(between the vector n and the velocity direction) 

0* = (l/2)arccos(l/D) . (35) 

Otherwise, for \D\ < 1, the tumbling regime takes place: 
the vector n experiences a time-periodic motion with an 
average rotation in the shear plane. Thus, D ~ 1 corre- 
sponds to the tank-treading to tumbling transition. As 
follows from Eq. ([33]), the transition is realized via the 
saddle-node bifurcation. 



IV. PERTURBATION EXPANSION: GENERAL 
DYNAMICS 

Here, we start to carry out the program (formulated 
in Subsection III C|) leading to a dynamic equation for 
the dimensionless displacement u. The program can be 
realized for nearly spherical vesicles by using a general- 
ization of the Lamb scheme. In accordance with Lamb 
[39I I (see also Ref. [i^l), a solution of the stationary 
Stokes equation can be explicitly expressed via the ve- 
locity field taken at a sphere both for the internal and 
external problems. The Lamb scheme can be directly 
applied to a spherical solid body immersed into a fluid 
or to a spherical cavity filled up by a fluid. Then the ve- 
locity field is expressed in terms of its surface value. For 
a nearly spherical vesicle the scheme has to be slightly 
modified. Namely, one should express the velocity field 
via its value on the sphere of the radius tq. The values 
can be obtained by an analytical continuation of the in- 
ternal and of the external velocity fields and are slightly 
different for the internal and for the external problems. 
The next idea is to represent the boundary conditions 
as an expansion over the displacement u which is small 
parameter for nearly spherical vesicles. 

In the zeroth approximation, one can ascribe the mem- 
brane velocity directly to the sphere r = tq ignoring devi- 
ations of the vesicle shape from the sphere. Keeping then 
the lowest in u terms in all expressions one passes to an 
equation for the dis plac ement u equivalent to the one dis- 
cussed in Refs. [2^ |25|. However, as we demonstrated in 
Ref. [4^, such approximation is not self-consistent. The 
problem is that it leads to a dynamics sensitive to initial 
conditions. And one can overcome this sensitivity only 
by accounting for high order terms in u. 

A. Closed equation 

Here, we derive an equation for the displacement u 
in the approximation where the membrane velocity and 
the boundary conditions (|13I14[) are related to the sphere 
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r = rp . Corrections to the equations associated with the 
deviations of the vesicle shape from the sphere are small 
in u. However, for the reasons formulated above, we keep 
the leading non-linear in u term in the expression for the 
boundary force Q. 

Note, first of all, that the variational derivative of the 
effective free energy (|2ip can be represented as 

— =ro{ -k[H{H^/2 - 2K) + A^H]+ aH] . 
du 

Therefore the boundary condition p3p can be rewritten 
as 

ST/Su = -25- + roPj„ - roPout- (36) 

Here, we divided the surface tension, a, into a homoge- 
neous, a, and an inhomogeneous, ct, parts. By definition, 
the zero angular harmonic is absent in a. Next, for the 
sphere r = tq the average curvature is H = 2/rQ and 
d^lnd^Vi (X S^d^Vi = 0, that explains validity of the 
expression (fSGj) . 

To find the inhomogeneous part of the surface tension, 
a, one has to use the second boundary condition, p4|) . 
Taking the derivative of the relation (fT4|) and relating 
a result to the sphere r = ?'o one obtains 

l{l + - 2ai + 2)il ~ l)vr,i 
= ij[{l + 2){l-l)vr,i+rld\j]^^ 

where cr; and u^,; are contributions to the surface tension 
and to the radial velocity, respectively, associated with 
the ^-th order angular harmonic. As above, the subscripts 
in and out are related to the interior and exterior regions 
of the vesicle. 

Applying the Lamb scheme to the sphere r = rg, one 
finds for the internal problem 

r' (?- l)(2Z + 3)^ , , 

/ '"o ' 

-r-}:[^^-^^)dtu, (39) 

Here, we used the condition drVr{ro) = 0, follow- 
ing from the relation UlkdiVk = where U = 
(sin0cos(y9, sin0sin(/3, COS0) is the unit vector perpen- 
dicular to the sphere r = r^. We substituted also 
Vr{To) = r^dtu, which is the kinematic relation (jl9p taken 
in the main approximation in u. 

For the external problem, one should separate the ex- 
ternal flow velocity V since it does not tend to zero as 
r oo. Writing v = V + w, one obtains from the Lamb 
scheme 

= L {-^i ^ - ^2 ) ^^^^ 



where, again, we used the incompressibility condition 
drWr{To) = and the kinematic relation Wr = r^dtu. 
However, the expressions (|40l4ip should be modified for 
/ = 2 because of the external flow. The modified incom- 
pressibility condition is drWr.2{i'o) + Sikhlk = and the 
modified kinematic condition is dtU2 — Wr.2 + foSikhlk- 
The modified conditions lead to the expressions 

Pout,2 = ?7(49tU(2) - 5s,kklk)rl/r^, (42) 
Wr^2 = (2i9tU(2) - 5sikklk/2) rl/r^ 

-{dtu^2)-SsM/2)ryr\ (43) 

Collecting together the relations (|36M3p one finds a 
closed equation for the displacement u 

1 SJ^ 

a{dt - u;d^)u = 10s,jhlj — . (44) 

Tyro Su 

where a is a dimensionless operator with angular compo- 
nents 

_ 2/3 + 2P + 3/2 - 5 t7 f + l-2 4C 

~ /(Z + 1) ^ /(/ + !) ~i ^ 1(1 + 1) 

We included into Eq. (gj]) a dependence on the rotational 
part of the external flow, which can be established by an 
account of the non-linear term in the kinematic relation 
(|19p . The Z-axis of our reference frame is implied to be 
directed opposite to the angular velocity a; of the external 
flow. 

Recall that the quantity a entering the dynamic equa- 
tion (UD) through Eq. (|2ip is the surface tension a aver- 
aged over angles. As previously, tr is an auxiliary quan- 
tity ensuring the surface conservation law. Let us stress 
that (T is a function of time adjusting to the current vesi- 
cle shape. Note that the strain and the rotation parts 
of the external flow are separated: the angular veloc- 
ity u) extends the time derivative (its effect is equivalent 
to passing to the rotating reference frame) whereas the 
strain matrix enters the term Sijlilj playing a role simi- 
lar to the free energy derivative. The reason is that the 
elongational part of the flow leads to some viscous dis- 
sipation whereas the solid rotation does not imply any 
dissipation. 

One can further elaborate the equation (P^ . First of 
all, one can keep in the effective free energy (PT|) terms 
of the second and of the third order in u. Higher order 
terms in T are small since u <^ 1. The reason why the 
third order term should be kept besides the second order 
one is explained above. Next, the expansion of the term 
Sijlilj does not contain any angular harmonics with I > 2, 
so this term does not push u outside the I = 2 subspace. 
Therefore the higher order angular harmonics in u die out 
after a flnite time due to the last term in ([44]) . Therefore 
one can use a reduced equation where u contains only 
second order angular harmonics. The operator a in this 
case is reduced to a constant 

16 / 23Ty C \ 
3 V 32 Ty 2Tyro/ 
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depending on the viscous ratios. The constant a can 
be called a generalized viscosity contrast. Note that the 
limit a ^ oo (where the viscosity of the internal fluid 
or the membrane viscosity tend to infinity) should corre- 
spond to a solid body behavior of the vesicle. 

One can substitute the expression (j26p into the equa- 
tion (|44p and then project the equation on the subspace 
1 = 2. Then the equation (|44|) is reduced to 



24 



(46) 







f 3k 


^ V 1 











where the subscripts in (9^ti)^ and {sijlilj)^ designate 
projections to the basis calculated in accordance 
with Eq. (|24)) . 



B. Rescaled equation 



The factor a in Eq. should be extracted from 

the condition 2ii^w^ = A, which is the area conservation 
law written in the main approximation in A. Then one 
obtains 



, , , , 24k „ 

l\){Sijlilj)p -I 3 ^p^xU^Ux 

Wo 



(47) 



where terms of the order su are neglected. In accordance 
with Eq. the right-hand side of Eq. (|Tf)) is a sum 

of two terms, proportional to Sij and to k. Note that 
the term proportional to k, is of the second order in u 
(the first order term is absent), that justifies keeping this 
high-order term in the expansion of the free energy. 

Deriving the equation (|T7)) we neglected terms of or- 
der su. They are much less compared to the term with 
K provided s ^ K-\/A/(??rQ). This is the applicability 
condition of the equation. However, below we demon- 
strate that results obtained from the equation (|17|) can 
be extended to stronger flows. 

After some rescaling the equation (f47|) can be rewritten 

as 



T, dt - ) Up = (Spp - UpUp)iSp + Ep,xUMx), (48) 



Up = V2up/VA, Ul + ... + Ul^ 1, (49) 

where S = (7-y/7r/A/5)S. The parameters in the equation 
are defined as follows 



5* = 



I47r srjTQ 



3\/3 kA 
where = SijSij/2. 



A = 



70r ar]rl 
12^/T0KVA' 
\/3 \fKauj 



(50) 
(51) 



The "vector" Sp in Eq. p5)) has an absolute value 
5, given by Eq. ([5T|) . its "direction" is determined by 
the projections of the object Sijklj to the basis ([23|) that 
is the "direction" is determined by the structure of the 
strain matrix s^ . Therefore the two parameters, S and 
A, together with the "direction" of the "vector" Sp com- 
pletely determine a character of the vesicle dynamics in 
the external flow. 

The quantity is the characteristic time scale of the 
vesicle relaxation. In comparison with the combination 
Tyrjj/K, related to the external fluid, the time r^,, contains 
an additional factor a, reflecting the contributions of the 
internal fluid viscosity and the membrane viscosity into 
the relaxation, and also the factor A^^/^. This extra 
factor reflects the slowness of the second order angular 
harmonic relaxation related to the degeneracy of the sec- 
ond order free energy in m. Due to the degeneracy the 
relaxation is determined by the third order term in u in 
the effective free energy, which contains a smalhiess A^/^ 
in comparison with the energy of higher angular harmon- 
ics. Therefore the adiabaticity condition, enabling one 
to use the stationary Stokes equation, can be written as 

^ QT^/rj that is arf'rQ ^ pn^/A. The inequality is 
valid because of the large value of the radius tq (in com- 
parison with the molecular length) and the small value 
of A. 

The parameter S characterizes the relative strength of 
the external flow. The expression ([5T|) for S can be ex- 
plained as follows. The external viscous surface force rjs 
should be balanced by the surface tension a times a varia- 
tion of the vesicle curvature which is estimated as ^/A/vq. 
Therefore a ^ rjsro / VA. As we have established in Sub- 
section [lir^ see Eq. (p9)) . the characteristic equilibrium 
surface tension can be estimated ctq ~ k-v/A/^q. Ratio 
of the surface tension values is therefore estimated as S: 
S ^ ct/cto- The applicability condition s <^ K,y/A/(rjr^) 
of the equations (|47l48p can be rewritten as 5* <C 1/ a/A, 
in terms of 5*. 

The parameter A determines the relative strength of 
the rotational part of the external flow. Note, that 
bJT^ ~ SK. The condition A ~ 1 corresponds to the 
angular velocity whose effect is comparable with the ef- 
fect of the strain, the condition gives uj ^ s/ {a^/A). This 
characteristic angular velocity does not coincide with the 
characteristic value of s, that stresses again different roles 
of the rotational and of the elongational parts of the ex- 
ternal flow. 



C. Very strong external flows 

Let us consider the case of very strong external flows, 
S > 1/VA, that is s > K^/A/(ryri^). As we already 
explained, the membrane surface tension, which is deter- 
mined by the balance between the viscous surface force 
jys and the product of surface tension a and the vari- 
ation of the vesicle curvature ^fK/rQ, is estimated as 
a ~ risro/y/A. Therefore for the external flows, we con- 
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sidcr here, a ^ i^/tq- The inequality implies that the 
leading role in the vesicle reaction to the external flow is 
played by the surface tension. 

Based on the inequality a 3> k/tq, one could try to 
neglect the terms with the module k in Eq. to ob- 
tain an equation independent of the curvature energy ((2|). 
That would correspond to omitting the contribution v^'^^ 
in the decomposition (jlSp . However, the resulting equa- 
tion is incorrect. First, there are additional terms in the 
equation for u, originating from the regular expansion 
over u and related to an account of deviations of the vesi- 
cle shape from the sphere r = tq. The additional terms, 
which can be estimated as su, are of the same order as 
kept in Eq. pS)) . We ignored the terms in the equation 



(|T7)) for u at s ^ Ky/A/{r/rQ) since then the term ^ su 
is small. Now the additional terms should be taken into 
account. Second, as we will see in the case of the planar 
flow, the truncated equation (without K-terms) cannot 
be used for describing the vesicle dynamics because of its 
internal properties. 

If the terms with the module k arc neglected in the 
boundary conditions (|13ll4p then they appear to be 
invariant under the transformation v —v, P 
—P, a — > —a. The stationary Stokes equation rjSI^v = 
VP and the boundary condition ([6|) are also invariant un- 
der the transformation. The kinematic relation (fTO|) be- 
comes invariant under the transformation if one adds the 
rule t —>■ —t. Therefore, in this approximation the back- 
ward in time evolution of the displacement u is equivalent 
to the direct evolution in the external flow with the veloc- 
ity —V. Since the vesicle dynamics is determined by the 
matrix ([5]), the transformation V —^ —V is equivalent to 
space inversion. That produces some additional symme- 
try leading to essential consequences for planar velocity 
fields. 



V. PLANAR VELOCITY FIELD 

In this section we discuss in more details the case of a 
planar external velocity field, when the velocity vector V 
lies in a plane and is independent of a coordinate normal 
to the plane. Then the velocity gradient matrix djVi has 
nonzero elements only in the plane being a 2d matrix. We 
choose X and Y axes of our reference system parallel to 
the plane and assume without the loss of generality that 
the diagonal elements of the matrix djVi are zero. Two 
non-diagonal components of the matrix dj Vi completely 
determine the flow, they can be parameterized in terms 
of the strain s and the angular velocity oj as dyVx = s-\-lij 



and dxVy 



Lo. In particular, for an external shear 



flow the only nonzero element of the matrix is dyVx ~ 7 
and LO = s = 7/2. 

For a planar velocity field, the only nonzero element of 
Sfj, in Eq. (|48)) is 5*5 = S. Thus only two parameters in 
Eq. S and A, completely determine the type of the 

vesicle dynamics. One of the main goals of this paper is 
to construct a "phase diagram" in the plane S — A, which 



indicates a type of the vesicle motion for a given pair of 
these parameters. 

Very strong flows, characterized by S* ^ I/a/A need 
a special consideration due to the reasons formulated 
above. It is presented in a separate subsection. Sur- 
prisingly, the case can be analyzed in the framework of 
the same equations (|47I48|) . 



General consideration 



For the planar flow, it is convenient to pass from the 



variables Un 



to another set of variables, "angles" 



<f>, 6, J, defined as 

Ui = sin e cos J, (52) 
U2 = sin O sin J cos 5*, U3 ^ sin 8 sin J sin 5", 
J74 = cos 8 cos(2$) , C/5 = cos 8 sin(2<I>) , 

where 8 varies from — 7r/2 to tt/2 and J varies from 
to 7t/2. The representation ([5^ satisfies the normaliza- 
tion condition (|49p and, correspondingly, consists of four 
variables instead of five components C/^ . 

As follows from the equations ([35]), there is a solution 
with J = 0. Actually, it is a consequence of general sym- 
metry 6 —9, which is present in the vesicle equations 
in the planar flow. Due to this symmetry, there exists a 
symmetric in 6 solution of u, that corresponds to J = 
as it follows from Eqs. (|23l25l52p . A special analysis is 
required to analyze the stability of the solution. Our nu- 
merical simulations demonstrated the stability. That is 
why we mainly limit ourselves to the investigation of this 
case. At J = the displacement u is 



cos 8 sin^ cos(2(^ - 2$) + sin 8 



1 - 3 cos2 e 

7! 



(53) 



Therefore the "angle" $ characterizes the vesicle orien- 
tation in the X — Y plane whereas the "angle" 8 deter- 
mines the vesicle shape. In the case J = the system of 
equations (gS]) is reduced to the following two equations 



S'sinesin(2$) + cos(38), 
cos(2$) 



cos 8 



A 



(54) 
(55) 



Note that the last, nonlinear in u, summand in (j47p pro- 
duces the term cos(38) in Eq. and no term in Eq. 

Let us flnd a region of parameters S and A where the 
equations ()54l55p admit stable stationary points, the case 
corresponds to the tank-treading vesicle motion. Equat- 
ing to zero the right hand sides of the equations, one 
flnds relations determining a stationary point (for given 
parameters S and A). In order to investigate its stabil- 
ity one should linearize the equations (|54l55p near the 
stationary point to obtain 



ndtiSej^) = B(J8,(5$). 



(56) 



10 



Tumbling 




damping relaxation 



10 



FIG. 3: Phase diagram 



and the trembling regime corresponds to a cycle separat- 
ing the atlas into two regions one of which does not con- 
tain any of the poles. In the tumbling regime the "angle" 
<i> grows without a limit whereas in the trembling regime 
it varies in a restricted domain. A transition line from 
tumbling to trembling, obtained numerically, is depicted 
by a dashed line in Fig. [31 

The diagram has a complicated structure near the spe- 
cial point S ~ a/3, a = 2/ a/3. A vicinity of the point is 
depicted on Fig. 21 where the regions of coexistence of 
two different stable points and of a stable point and of a 
limit cycle are shown. More detailed description is given 
below. 

To avoid a misunderstanding, note that there is an 
additional region in the S — K plane where stationary 
solutions exist, which are stable in terms of the vari- 
ables and $. However, a stability investigation in 
the framework of the complete equation ((48|) shows that 
these solutions arc unstable in the extended space with 
the variables J and Therefore these solutions cannot 
be realized as the tank-treading motion in real systems. 




El 



1.72 1.74 1.76 1.78 1.8 1.82 1.84 1.86 1.88 1.9 

S 

FIG. 4: Phase diagram in a vicinity of the special point. 



The stationary point is stable, if both eigenvalues of the 
matrix B have negative real parts. Thus the stability 
conditions are tr i? < and det S > 0. A region in the 
S-K plane where stationary stable points exist is shown 
on the Fig. [31 The region below the line A = 2/-\/3 
corresponds to positive values of <&. The tank-treading 
regime above the line A = 2/\/3 leads to negative values 
of $ for 5* > V3. 

For the upper region in Fig. Othe attractors of the sys- 
tem (|54l55p are limit cycles. They correspond to either 
tumbling or trembling behavior. The difference is illus- 
trated in Fig. O where the — $ atlas is plotted. The 
tumbling regime corresponds to a cycle separating the 
atlas into two regions each containing a pole 8 = ±7r/2. 



B. Tank-treading to trembling transition 

The tank-treading to trembling transition is deter- 
mined by the condition iv B — 0. The corresponding 
curve on the S — K plane starts from the above special 
point S = V3, A = 2/-\/3 (the point ei on Fig. |4]) and 
goes to the right. This curve, marked as red, is described 
by the equation 

A* = V2v/1 - l/S'2, (57) 

where S varies from \/3 to oo. Above the red curve, 
at A > A*, the stationary point loses its stability via a 
Hopf bifurcation. Let us establish characteristics of the 
bifurcation. 

Expanding the equations (|54l55p near the point (|57p . 
one finds the equation for a complex variable z 

ndtz = ez- iy/S^ - 3 z - K\zfz, (58) 

where 

8SVS^.. ^ 2V2VS^^{S^ + 5) 
^-^^^(A-A), K= , 3)3/2 • 

The variable z is expressed via the deviations of the "an- 
gles" from their stationary values as 



Above the transition line, at A > A*, the vesicle mo- 
tion is described by a limit cycle with the radius propor- 
tional to V-A — A*, near the transition curve. This mo- 
tion corresponds to trembling since the radius is small. 
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and the corresponding limit cycle cannot surround a pole, 
see the atlas in Fig. [S] 

Note the critical dependence of all parameters in Eq. 
([58|) on 5 — Vs. Taking into account the critical depen- 
dence, one concludes that near the line the amplitudes 
of the 8 and $ variations can be estimated as ^/A — A*, 
without a critical dependence on S" — v^- 

A vicinity of the special point ei needs an additional 
analysis since the frequency t~^\/S^ — 3 of the Hopf bi- 
furcation tends to zero at the point and the approxima- 
tion leading to the equation ([55)1 is not valid there. 



C. Tank-treading to tumbling transition 

The transition from tank-treading to tumbling is de- 
termined by the condition det B — 0. The transition 
curve on the S ~ A plane, designated as orange, has a 
complicated shape, it can be described in a parametric 
form 



for the parameter ^ which is 

T,dt^^5AFi+F2^'', 



(60) 



S 



CVl5-32C2-^-16C4 



A = 



(59) 



where the parameter C varies from 1/V2 upto V3/2. The 
boundary value C = 1/ corresponds to the above spe- 
cial point S ~ VS, A = 2/\/3, and the boundary value 
C = corresponds to the point S = 0, A = l/\/3. 

The expression for S has a maximum at Co = Vl — 2^'^/^, 
for C close to Co we obtain S « 1.8737 - 46.97(C - Co)^- 
The value C = Co corresponds to the turning point cq on 
Fig. HI 

To be more precise the condition det B ~ determines 
a stability boundary of the tank-treading regime. The re- 
gion of parameters l/\/2 < C < Coi determining the part 
of the orange curve going from the point ei to the point 
eo, corresponds to a transition from one tank-treading 
regime to another one. That is why the dark green region 
on Fig. m corresponds to coexisting two tank-treading 
regimes. At passing from the region through the red 
curve the Hopf bifurcation occurs. However, the bifur- 
cation takes place for one of two possible tank-treading 
regimes, the other one remains stable. Therefore there 
exists a region of coexistence of the tank-treading and 
trembling, the region is marked as fuchsia on Fig. U) 
Its left boundary, corresponding to an instability of the 
trembling motion, is found numerically. 

Let us consider now the region of parameters Co < C < 
•\/3/2, giving the upper part of the orange curve. We are 
interested in the dynamics of the deviations (56. (5$ from 
the stationary values of the "angles" . An analysis shows 
that there exists a linear combination C of the deviations 
possessing a slow dynamics, that is the characteristic re- 
laxation time of C scales as VSA where i5A = A — A{S) 
and A(S') is the value of the parameter A at the transition 
curve, determined by the expressions (j59p . Then, using 
adiabaticity, it is possible to formulate a closed equation 



valid at small 6 A and ^. The expression ([50]) is charac- 
teristic of a saddle-node bifurcation. 

The parameters of the saddle-node bifurcation ((60|) 
have critical behavior near the boundary points C = Co 
and C = a/3/2 which can be expressed as 



i^i =/iv/3^4?, = /2V3^4C^(C"Co), (61) 

where /i and /2 are functions of C varying less than by 
25% as C runs from Co to V3/2. In the normalization 
where ^ (^e)2 + (j$)2 g^^^g ^ 4 5 ^^^^ ^ qq 

at C ^ Co- The function F2 tends to zero as C ^ Co- 
Therefore higher order terms in the equation for C should 
be taken into account near the point. 

The limit cycle which is a result of the saddle-node bi- 
furcation, destroying the tank-treading regime, could be 
unstable in its turn. The situation is realized between 
the points eo and 65, see Fig. |4l Above the segment 6064 
a final result of the instability is another tank-treading 
regime continuously continuing to larger S. Above the 
segment 6465 it is trembling, characterized by a limit cy- 
cle which does not pass through the stationary point. 
After the point 65 (to the left from the point) the tank- 
treading regime is destroyed and a limit cycle passing 
through the stationary point is formed. Above the seg- 
ment 6560 the cycle corresponds to trembling, otherwise 
it corresponds to tumbling. 

Both functions, Fi and F2, tend to zero as C ^ \/3/2. 
This is the limit of weak external flows, where the vesi- 
cle relaxation rate is proportional to the strength of the 
flow. Then the right-hand side of the equation should 
be proportional to S which behaves like S oc •\/3 — 4C^ 
as a consequence of Eq. ((59|) . That explains the depen- 



dence Fi,F2 cx ^3 — 4C^. However, the limit of weak 
flows needs an additional analysis since, in accordance 
with results of Subsection llll B[ there are two soft degrees 
of freedom corresponding to rotations of the equilibrium 
uniaxial ellipsoid. We postpone the analysis to the next 
section. 



VI. LIMIT CASES 

We established a general picture of the vesicle dynam- 
ics in an external flow which appears to be rich of dif- 
ferent types of behavior. Particularly, the phase diagram 
depicted in Fig. [3] contains different domains and has 
a complicated structure. The situation is simplified for 
different limiting cases which can be analyzed in more 
detail. 

For example, the equation ([55]) leads to the conclu- 
sion that in the case A ^ 00 the vesicle rotates with the 
angular velocity uj. For external flows with comparable 
strain and rotational parts vesicle rotation is quite natu- 
ral, since in accordance with the definition (|5ip the limit 
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A ^ oo is achieved either at a ^ oo or at s — > 0. The first 
case corresponds to a sohd body behavior of the vesicle, 
so one should reproduce the classical Jeffery's result [46j . 
which predicts, that for external flows with uj > -v/As 
the tumbling regime supersedes the tank-treading one. 
The second case corresponds to a purely rotational ex- 
ternal flow, where the fluid rotates as a whole with all 
inclusions. 

Below, we analyze more complicated limit cases. 



A. Purely elongational flow 



The purely elongational flow is realized provided the 
angular velocity oj is equal to zero, that is dyVx = dxVy = 
s. Therefore, in our designations, the elongation is di- 
rected along the main diagonal in the X — Y plane. 

The condition lu ~ Q leads to A = 0, in accordance with 
the definition ([5T|) . In this case the system of equations 
(|54l55p has a stable stationary point (f>o, Go, determined 
by the relations 



$0 = 7r/4, 5'sin9o = cos(39o). 



(62) 



The "angle" 8o monotonically decreases from 7r/6 to zero 
as S increases from zero to infinity. The value $o = 
7r/4 is quite natural since it corresponds to the vesicle 
orientation along the elongation direction, as is seen from 
Eq. dSSl). 

The stability of the point ([5^ can be easily established 
from the linearized equations 



cos Ho 

ndtSe - [S" cos Go + 3sin(3eo)] 50. 



(63) 
(64) 



In the limit 5* ^ 1 both "angles" relax to their equilib- 
rium values with the same rate SVIOtt s(3'\/3A a)^^. 

Recently a wrinkling effect was observed in purely elon- 
gational flows at a sudden invertion of the elongation di- 
rection, see Ref. [i^. The effect can be explained in the 
framework of our scheme, the corresponding analysis is 
presented in Ref. [sot . 



B. Weak external flows 

Let us consider weak external flows characterized by 
the condition 5* <C 1. We have already discussed the 
case in Section IIII Bl from the phenomcnological point of 
view. We can analyze the case in terms of the "angles" Q 
and $ and then establish a value of the phenomenological 
constant D, introduced in Eq. (PT|) . 

As follows from Eq. (O, for S* < 1 the "angle" 6 
is close to 7r/6, which is a stable point of the equation. 
Substituting the value O = 7r/6 into Eq. ([55]) one obtains 
a closed equation for the "angle" $ 



r,at$ = (5/V3) cos(2$) - S'A/2. 



If A < 2/\/3 then the equation ((65)) has a stationary 
point 



$ = - arccos I — - — 



(66) 



which is stable. Otherwise $ increases unlimited, that 
corresponds to the tumbling regime. Therefore A = 2/V3 
is the transition point from tank-treading to tumbling. 

If 6 = 7r/6 then the expression ([55]) describes a prolate 
uniaxial ellipsoid with the principal axis directed along 
the vector ([32)) with (/> = $ and = 0. Comparing then 
the equation ([55)) with Eq. ([55)) (obtained for a shear 
flow with s = Lo = 7/2) one flnds 



D = 



3aVA 



(67) 



As it should be, the transition point D = 1 from tank- 
treading to tumbhng corresponds to A = 2/ a/3. Let us 
stress that the value ([67)) is independent of the character 
of the external flow. Therefore the equation ([5T)) with 
([S7)) is correct for any external flow V. 

Note that in the "solid body" limit a ^ 00 (where the 
viscosity of the internal fluid or the membrane viscosity 
tend to infinity) the quantity D tends to zero. Account 
of higher order terms in A gives, that in the solid limit 
D stops decrease at value order of -s/A. Diminishing of 
D leads to a solid rotation of the vesicle in particular 
case of external shear flow as follows from Eqs. ([33l34p . 
The behavior corresponds to the classical result of Jeffery 
[iHl, who demonstrated that a solid ellipsoid rotates in 
an external planar flow, provided lu > \/As. 



C. Strong external flows, truncated equations 

In the case of strong external flows, where S is large, 
the last (non-linear in U) term in the right-hand side of 
Eq. ()48p is small in comparison with the first one. If 
to omit this last term we pass to a truncated equation. 
In terms of the variables introduced by Eq. ([5^ the 
truncated equation is written as a system of equations 

(n/S')(9te = -sinesin(2$). 



(n/5)9f$ = 



1 



cos(2$) 



cos 8 



-A 



{T,/S)dt^^ A/2, 



dtJ = 0, 



(68) 



(65) 



homogeneous in S. The system (1681) corresponds to the 
limit case considered by Misbah [2J] and Vlahovska and 
Gracia [2^ . In the subsection we examine solutions of 
the system ([68| . A relation to observable behavior of 
vesicles, which is not straightforward, is discussed in the 
next subsection. 
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FIG. 5: Vesicle dynamics on the O — $ atlas. 

The system (|68p leads to conservation of two quanti- 
ties, J and an additional integral T, which can be intro- 
duced via the relation 



A-cosT A-cosecos(2$)' ^ ' 

For definiteness, we choose a root of the equation 
lying in the domain |T| < arccos(l/A). Existence of two 
integrals of motion, J and T, implies that a character 
of an evolution, described by the system ([68]) . depends 
on initial conditions (determining the values of the inte- 
grals). 

Thanks to existence of two integrals of motion, the 
system of equations (|68)) can be completely integrated. 
For the purpose we introduce a variable 

p = exp js/r, ^ At' t/5(i')| • (70) 

It turns out, that for arbitrary initial conditions a solu- 
tion passes through a point, where = 0. It is conve- 
nient to choose initial time, t = 0, as a moment, when 
?75 = and [/4 = cosT > 0. Then the initial conditions 
are p = 1 and dtp ~ Q and one derives from the system 
the following equation 

(r,/5)2 9,V= -(A'-1)p + A2 - AcosT, (71) 

which can be obviously solved explicitly. The parameters 
$ and O are expressed via the variables p and T as 

[/4 cos e cos(2$) = Ap - A -I- cos T (72) 
C/s = cos e sin(2$) = {T^/S)dtp/p. (73) 

A solution of the equation ([7T|) behaves differently at 
A < 1 and at A > 1. At A < 1 the variable p trends 
to infinity as time grows, thus C/5 has its utmost value 
be equal to Vl — A^, that corresponds to a stable point 
(tank-treading behavior). At A > 1 the variable p experi- 
ences oscillations. Then one conclude from Eqs. (|71H73I) . 
that the vesicle evolution is described by a limit cycle, its 



characteristics are determined by values of the integrals 
J and T. If the angle $ grows (decreases) unhmited, the 
vesicle is in the tumbling regime. On the contrary, if <f> 
is bounded, one deals with the trembling regime. It is 
convenient to represent vesicle dynamics as a geographic 
atlas, where 9 and $ play roles of latitude and longi- 
tude, correspondingly, see Fig. [51 The trembling regime 
corresponds to a closed curve on the Q-^ atlas, which 
does not surround a pole. The tumbling regime corre- 
sponds to a curve separating the poles. In Fig. [5] such 
curves terminate at the boundaries of the atlas having 
the same Q on the ends, since points on the right and on 
the left boundaries of the atlas with the same latitude are 
physically identical. For the truncated system ([68]) the 
tumbling and trembling regimes coexist at any A > 1. A 
choice between the regimes is determined by a value of 
T. If cosT > 2A/(A2 + 1), then the limit cycle corre- 
sponds to tumbling, otherwise the cycle corresponds to 
trembling. If T takes one of its boundary values, that is 
if cosT = 1/A, then the limit cycle degenerates into a 
point 

$ = 0, cose = 1/A. (74) 

Thus, the truncated system of equations (|68|) has two 
stationary points corresponding to the tank-treading 
regime. 

D. Strong external flows, slow dynamics 

We demonstrated that the truncated system of equa- 
tions ([68)) can be completely integrated. However, the 
system cannot be directly used for an analysis of the vesi- 
cle dynamics in the limit of strong shears 5 3> 1. Indeed, 
the system leads to a dependence of its solution on ini- 
tial conditions and admits different limit cycles for any 
A > 1 . Both these properties contradict obviously to the 
results obtained in Subsec. IV Al The contradiction is re- 
solved if one restores the terms omitted in the truncated 
system (|68p and originating from the last (non-linear in 
U) term in the right-hand side of Eq. pS)) . The restored 
terms provide a relatively slow evolution of both inte- 
grals of motion, J and T, which leads to a well defined 
behavior, independent of the initial conditions. 

Our nearest goal is to deduce the equations of motion 
for J and T. We consider the case A > 1 where the 
truncated system of equations ([68]) leads to limit cycles. 
Then, analyzing the complete system of equations, one 
can separate fast motion along the limit cycles and rela- 
tively slow evolution of the integrals of motion on scales 
larger than the cycle period. Note, that a typical time of 
the fast dynamics is t^/S whereas a typical time of the 
slow dynamics is r*, the large ratio of the times justifies 
the separation. The equation controlling the slow evolu- 
tion can be found by averaging over the cycle period of 
the expressions for the time derivatives of J and T ob- 
tained from the complete system ([JSj). At this averaging, 
one can use the fast dynamics described by the system 
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([68)1 . The result can be schematically represented in the 
form dtT = T(T, J,A), dtJ = j(T, J,A). Explicit ex- 
pressions for T and J are quite cumbersome, so we do 
not present here its final form. 

One can check, that the ultimate value of J is equal to 
zero at any A. That is why below we consider the case 
J = 0, and find the slow dynamics for the quantity T 
at the condition. Then T is a function of A and T, the 
fmiction can be found from Eqs. (|54I55I69|) 

T = ((aT/ae)cos(3e)) , (75) 

where angular brackets mean averaging over the cycle 
period. The evolution of T described by the equation 
dfT = T leads to a stationary value which can be found 
from the equation T = 0. The value corresponds to 
a limit cycle, which is stable provided dT/dT < 0. 
In accordance with the analysis, made in the previous 
subsection, a character of the limit cycle depends on 
the value of T. If cosT > 2A/{A'^ + 1), then the 
limit cycle corresponds to tumbling, otherwise the cy- 
cle corresponds to trembling. A numerical investigation 
based on Eq. ([75]) shows that the boundary value of T, 
T = arccos[2A/(A^ -I- 1)], is achieved at A sa 1.52, larger 
A correspond to tumbling, smaller A correspond to trem- 
bling. 

At A < all the limit cycles appear to be unstable. 
That means that the condition the stationary point ([74|) 
appears to be an attractor. Thus, the system falls to the 
stationary point that corresponds to the tank-treading 
regime. Note that the position of the stationary point 
is slightly shifted due to the presence of the additional 
term in the equation for Q The new position is 

e = arccos(l/A), * - (76) 

Zij Sm \y 

in the main approximation in 1/S. A correction to the 
value O = arccos(l/A) is of the order l/S*^. 

E. Extremely strong flows 

The above analysis is, strictly speaking, correct for 
flows characterized by S" < 1/y/A. For stronger flows, 
1/\/A < S, our consideration should be extended. Some 
additional terms of the higher order in u should be taken 
into account in the equation for u, which for l/\/A < S 
are larger than those kept in Eq. ((47|) . Leading terms 
of such kind can be estimated as su, they originate, say, 
from account of deviations of the vesicle shape from a 
spherical one. The terms are associated with the contri- 
bution v^"^ to the membrane velocity, whereas the second 
term in the right hand side of Eq. ((47|) is associated with 
the contribution d^**' to the membrane velocity, see Eq. 

However, there is an essential difference between the 
terms t?^*-* and v^'^\ If one omits the term v^'^'^ then there 
appears a symmetry of the vesicle dynamics expressed 



in terms of the equation for u, which is invariant under 
simultaneous time and space invertions, see Subsection 
IIVCI Since the axes of our reference system are attached 
to the eigen vectors of the strain matrix s, the space 
inversion is equivalent to the transformation tp — *■ —ip, 
9^9, that is it can be written as $ ^ — $ and ^ — 
in terms of the "angles" ([52|) . Therefore the equations for 
Q and $ should be invariant under the transformation 

t^-t, $ ^ e ^ e. (77) 

One can easily check that the truncated equations ([68|) 
are invariant under the transformation (|77|) . However, 
our analysis demonstrated that the symmetry survives 
even though higher in u terms will be taken into account 
provided k ^ 0. 

If small in u corrections to the truncated equation (|68p 
will be taken into account then the system of limit cy- 
cles, characteristic of this equation at A > 1, will survive 
with small perturbations. Next, if the term v'^'^^ in the 
expression (fTS)) is neglected all the limit cycles will be 
neutral (they are no stable, no unstable). Indeed, the 
cycle passing through a point <& = 0, O cannot be stable, 
since then due to the symmetry (j77p it should remain sta- 
ble after the time inversion. Thus, the corrections to the 
truncated equation (|68p keep its main property which is 
an existence of an additional integral of motion (leading 
to the neutrality of the limit cycles). 

Thus only the term t;''*^-' proportional to n can destroy 
the integrability corresponding to the neutral limit cy- 
cles. It produces a selection leading to stable limit cy- 
cles or to stable stationary points. Since the selection is 
produced among limit cycles slightly disturbed in com- 
parison with the ones corresponding to the truncated 
equation ([55]) . the results will be the same in the main 
approximation in A. Thus the above results obtained 
for strong flows can be immediately extended to the ex- 
tremely strong flows, where S > 1/\/A. 

VII. CONCLUSION 

We have investigated dynamics of nearly spherical vesi- 
cles in an external stationary flow, being mainly focused 
on the shear flow. The general calculational scheme de- 
veloped for nearly spherical vesicles enabled us to analyze 
the dynamics in detail. The scheme is based on solving 
the id hydrodynamic (Stokes) equations with boundary 
conditions posed on the membrane. Besides the mem- 
brane bending elasticity we have taken into account the 
internal membrane viscosity, leading to an additional dis- 
sipation mechanism. 

There are essentially different regimes of the vesicle 
dynamics dependent on a relation between the strain of 
the external flow and the vesicle relaxation rate. In the 
weak flows the vesicle shape is close to an equilibrium one 
which is a prolate ellipsoid. Then a role of the external 
flow is reduced mainly to an orientation of the ellipsoid. 
In the strong flows the vesicle shape and orientation are 
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determined by the flow. We established the vesicle behav- 
ior for different strengthes of the external flow including 
its orientation relative to the external velocity. 

The most interesting phenomenon in the vesicle dy- 
namics is transition to tumbling regime which occurs at 
increasing the generalized viscous contrast ([^S]) or at in- 
creasing the rotational component of the external flow. 
We demonstrated that in weak flows a direct transition 
from tank-treading to tumbling occurs whereas in strong 
flows an intermediate regime, trembling, is realized. The 
behavior is in accordance with experiment of Kantsler 
and Steinberg Q and corresponds qualitatively to numer- 
ics of Noguchi and Gompper [2O] . The phase diagram we 
obtained is plotted in Fig. [3] where the variables S and 
A are some combinations of observable quantities (|5ip . 

The possibility of an intermediate regime between 
tank-treading and tumbling was discussed theoretically, 
first by Misba h [241 and then (qualitatively) by Noguchi 
and Gompper [20] ■ Note, however, that the calculational 
scheme used by Misbah [2J] and then by Vlahovska and 
Gracia [2^ is not self-consistent though formally the au- 
thors have taken into account principal terms of the equa- 
tion for the vesicle distortions. The situation cannot be 
improved by introducing higher order terms related to 
the external flow. And only after introducing higher or- 
der terms related to the membrane bending elasticity into 
the equation it enables one to find, say, boundaries of the 
region where the trembling regime is realized. 

The transitions tank-treading to tumbling and tank- 
treading to trembling have essentially different scenarios. 
The first transition is described as a saddle-node bifurca- 
tion whereas the second one is described as a Hopf bifur- 
cation. Our theory predicts existence of a special point on 
the diagram where the two above transition lines merge. 
One expects an essential "critical" slowness of the vesicle 
dynamics near the point. Therefore it plays a role anal- 
ogous to some extent to the critical point on fluid phase 
diagram. 

Near the transitions thermal fluctuations are relevant 
which smear the transitions. One expects that the effect 
is especially strong near the special point on the phase 
diagram. Role of the thermal fluctuations, which can be 
examined in the spirit of the works [i^, O, , consti- 
tutes a subject of special investigation, to be done sep- 
arately. Here we note only, that due to fluctuations one 
should be careful in comparison an experiment with the 
phase diagram obtained in our work since the diagram is 
deduced ignoring the fluctuations. 

It is worth to mention the effect, related to the ther- 
mal fluctuations, recently discovered by Kantsler et. al. 
p9| . It was shown that the relaxational dynamics of a 
vesicle in external clongational flow is accompanied by 
the formation of wrinkles on a membrane. Theoretical 
investigation of this effect presented in [HO] was based 
on the theory, developed in this paper. It was shown 
that the formation of wrinkles is related to the dynami- 
cal instability induced by negative surface tension of the 
membrane. 



We investigated nearly spherical vesicles assuming that 
the excess area factor A is small. It enables one to formu- 
late a powerful calculational scheme enabling to find de- 
tails of the vesicle dynamics analytically. We have doubts 
that the scheme can be generalized for a general case 
A ~ 1. Probably, the case can be investigated only nu- 
merically. However, the general approach should be the 
same: one has to solve the stationary Stokes equation 
inside and outside the vesicle with the boundary condi- 
tions (|6ll3ll4p at a given vesicle shape and then to use 
the equation (|19p to formulate an equation for membrane 
distortions. 
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APPENDIX A: GENERAL DYNAMICS 

In the appendix we present the derivation of Equation 
(|44p . For brevity we introduce viscosity contrast param- 
eter \ ^ fj/rj and surface viscosity parameter /i = C^/rjrQ 

At low Reynolds numbers velocity field of fluid is com- 
pletely determined by the boundary conditions. One can 
solve the Stokes equation @ inside and outside the vesi- 
cle membrane with boundary conditions on the mem- 
brane and a fixed boundary condition at the infinity ([5]) . 
This allows one to exclude the fluid motion and arrive 
to the close equation onto membrane surface. On the 
way one obtains integral equation, describing membrane 
surface motion [l3, Ej. Further analytical progress is 
possible only at the limit A ^ 1, when it is possible to 
expand the equation describing the membrane dynamics 
in series over A. This expansion has no regular limit at 
A ^ 0. Physically, exactly spherical vesicle can not con- 
serve its volume and area in arbitrary small external flow: 
to describe the vesicle behavior in an external flow one 
should restore final surface compressibility of the vesicle 
membrane. As a consequence, expansion series of some 
quantities, which depend on the vesicle shape, starts from 
zero or even negative order in A terms. 

First derivation of this type of equation was presented 
in |22| in particular case of equal inner and outer viscosi- 
ties, 77 = 77 and zero surface viscosity, C = (see also 
further generalizations [l^, HH], where case A 7^ 1 was 
considered). To proceed to the derivation, let us first 
exclude fluid motion inside and outside the vesicle. We 



parameterize velocity field inside and perturbed part of 
velocity field v — V outside the vesicle by the functions 



(Al) 
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correspondingly. The functions (jAl[) are functions of 
spherical angles 9, ip. Procedure of restoration of ve- 
locity field into the bulk from the functions is known as 
Lamb solution. The physical meaning of the functions 
(|Aip is the following. On a sphere with the radius the 
components of velocity and vorticity fields are given by 



r^drv' = r'", row'' 



(A2) 



where a; is vorticity. Using the Lamb solution one can 
make analytical continuation of the velocity into whole 
interior of the vesicle. We expand each of functions (|A1[) 
into series over spherical harmonics: for example, 



(A3) 



i>0,m 



where m runs from — / to /. Assuming that the velocity 
field is a regular function inside the sphere, one finds by 
solving of (III) that 



V = gradll + rot[rn7i] + rot rot [r-n72] 



where 



{I + 1)X'"'' - Y' 



(ro) 



/-I 



W2 



tUi 



E 



E 
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{roY 1(1 + 1) 



(ro) 



2l{l + l) 



(A4) 



(A5) 



(A6) 



(A7) 



In (jA5IIA7P the summation goes over m as well, although 
we omit it in the equatoins for brevity. The pressure field 



r' {21 + 3) (r'"'' -{I- i)a:'"^') 



I 



is induced only by the last velocity component in (|A4p 



Also we introduce covariant derivative on the surface V^, 
and for any quantity, for example a, we use the notation 

Force balance equations p3ll4p . written in the curvi- 
linear coordinates, have the following form: 



'a/3 



Ha-K [H{H^/2 - 2K) + A^i/] = 



(A9) 



where S is doubled symmetric part of velocity gradient, in 
Cartesian coordinates = diVj+djVi and brackets [. . .J 
denote the difference between the "inside" and "outside" 
values. The index "m" near vertical line indicates, that 
the equations are written on the membrane surface. 

We choose the two spherical angles as the internal coor- 
dinates, see (fT6)) . While developing a perturbation the- 
ory over A, it is convenient to use the method of do- 
main perturbations, and to use the sphere with the radius 
To, parameterized by the spherical angles, as intermedi- 
ate surface. Each point on the membrane surface corre- 
sponds to a point on the sphere with radius tq, having 
the same spherical angles. To distinguish vector quanti- 
ties and equations, defined on the intermediate surfaces, 
we write index "s" near vertical curve instead of index 
"m" for the membrane surfaces, see for example (|A9p . 
Any quantity defined in a bulk or on a membrane sur- 
face can be represented as a taylor expansion through 
its value and its derivatives defined on a intermediate 
sphere. This method is well defined for quasi-spherical 
membranes with A ^ 1. 

For the sphere we introduce one more object - direc- 
tional derivative, defined by the equation 



9a(} 



(AlO) 



where e^'^ is antisymmetric unit symbol, e^^ = 1. Also 
we denote Beltrami-Laplace operator on the sphere as 
O/rQ, such that O is Beltrami-Laplace operator on unit 
sphere. 



Curvilinear coordinates 



Utilization of Lamb solution. 



In the Appendix we use the curvilinear coordinates 
for the vector quantities in contrast to the main text, 
where we used cartesian coordinates for clarity. Let us 
parametrize the vesicle surface by two internal coordi- 
nates Leaving indices i, j, fc, . . . for vector quantities, 
written in cartesian coordinates, we identify vector quan- 
tities, projected onto the membrane surface and written 
in the internal coordinates by indices a,f3, . . .. We denote 
the metric and curvature tensors as gaf3 and hap corre- 
spondingly. Recall, that by definition the mean curva- 
ture of the membrane is given hy H = —g°'^hafj, where 
as usual we assume summation over repeating indices. 



In this subsection we express a list of quantities, which 
are involved in process of satisfying of boundary condi- 
tions, through functions (|Aip . We use the following nota- 
tions. Suppose fimction 17 is a function of spherical angles 
9, If and K(/) is some function of angular harmonics /. 
Then linear operator is defined as 



K+5 - E K(0 



9L 



(All) 



l>l,m 



Also we introduce linear operator K , whose definition 
can be obtained from (jAlip through replacing / ^ — (/ -I- 
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1). Also we denote 

K = K+-K", Q = AK^ 



(A12) 



We need in pressure, velocity and some components of 
symmetric part of velocity gradient Sy on sphere with 
radius rg. Here we list all quantities on inner side of the 
sphere, omitting index in. Quantities on outer side of 
the sphere can be obtaining by changing I ^ —I — 1 and 
plugging functions (jAip with index out. Pressure inside 
is 



p = roz/[tP^+ X + KPy+ Y] 



(A13) 



where 



il) = J^-'f^'\ ^ (A14) 

The velocity is given by 

Va = -Wid-^{2X + Y)-Wld-^Z. (A15) 
Tangential-radial components of velocity strain § is 



- V! 



■rsZ-\- 



In (jA16p necessary for us coefficients K are 



Tangential-tangential components of § are 
Sa0 = Uxg^p - [a-i [AX + 2Y] 



:q/3 



(A16) 



(A17) 



In the interesting case oi Y = Z = one obtains from 



= -2Vt 


1 + 20-1 


X 


S 







and 



/i^X^sL = 0. 
The radial derivative is given by 



(A18) 



(A19) 



(A20) 



where we kept the only necessary contribution from the 
unction X with 



fSX 

rrr 



(0 = 2(1-P). 



(A21) 



3. Obtaining of the evolution equation in main 
approximation in A. 

In this Subsection we return to the derivation of Eq. 
We have seven unknown scalar functions, six in 
(|Aip and surface tension a, which depend on spherical 
angles 6*, (p. To find these functions, one should satisfy 
all the boundary conditions on the membrane. Three 
boundary conditions come from velocity continuity in 
the whole space. For us it is convenient to use three 
equivalent conditions: continuity of normal component 
of velocity Ti;*, continuity of normal derivative of normal 
component of velocity n^dkv'' and normal component of 
vorticity I'^uj'' . Another three conditions come from the 
continuity of the momentum flux, (jl3ll4p . Seventh con- 
dition corresponds to the surface flow incompressibility, 
sec (jni). After the seven scalar filed are found, one should 
use the relation between the velocity and the temporal 
derivative of function u (fTB|) presented at p9p and ob- 
tain required equation of motion. In this section we put 
K = 0, and do not take into account the terms arising 
from the bending force of the membrane. 

To obtain the dynamical equation with the accuracy 
up to n-th order of v^, one should satisfy all boundary 
condition with the same accuracy. It is convenient to 
find the next in -v/A correction to the equation using the 
recursive procedure. In accordance with the scheme we 
represent any quantity, for example X'", as a series 



X' 



n=0 



Xl 



(A22) 



where oc A"/^. For vector quantities which have 
low indices we put the index n on the left from the main 
letter to avoid mixing of different indices. 

We represent the local value of the surface tension as 
a ~ a + a, where a does not depend on angles 6*, tp and 
the integral of a over the angles is zero. The reason of 
the division is different scaling laws of parts a and a with 
A: expansion of the quantities in series over A are 



(A23) 



)i=-i 



ri=0 



where lower index n near a term corresponds to the term 
scaling as A"/^. 

Let us find equation of motion in zero order in VA- 
Surface incompressibility, i.e. continuity of V'v^ reads 



where X'^ 



h, definition of h is 



Continuity of Pl^dkV^ leads to 



J n 



■yout ^ ycxt 



(A24) 



(A25) 



(A26) 
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where again y®'^* = h. For continuity of normal compo- 
nent of vorticity one should assume 

Z^" = + Z"^*, (A27) 

where Z'^^*' = —rQC^^^nPc,^^ and n is unit vector on sphere 
having radius tq. Membrane incompressibility condition 
yields the equation 



n 



0. 



(A28) 



Relation between the velocity field and the dynamics 
of the vesicle shape has the form 



X- = u. 



(A29) 



In order to derive the dynamic boundary conditions 
(fT3|l and (fT4|) . we should use the expansion of local mean 
H and Gaussian K curvatures in ^fK: 



2 



Hi = - 



{2 + d)u 



Ho = 



'"0 



1 



(A30) 



(A31) 



Now we have determined all the objects required to 
write down (|13p and (jl4p in the main approximation over 
A 



<j-iHi+aoHo=lPo\, (A32) 
+ Vtao= L^yoSraJI . (A33) 



Values of first terms in (jA32[) and (|A33p should be taken 
from (jA19p and (jAlSp correspondingly. In the approxi- 
mation Zq"' ~ 0, it is convenient to obtain the statement 
using the directional divergence (jAlOp of Eq. (|A33p : the 
resulting equation is a linear homogeneous equation on 
Zg"*. Taking simple divergence of Eq. (|A33p . one can 
find alternating part of the surface tension 

1 . . 5 



— 2 ^0 
Vro 



h. (A34) 



{2m(1 + 20-1)+Q-}w 2" 
Substituting (|A34p it into (|A32p one eventually obtains: 



uo 



Wo 



a(/) = 



_ 2(A+l)i3+3(A+l)i 



4/i (1 4 

5A+4 



20-1) , 

^ ^z-' /(i+l) ■ 



(A35) 



(A36) 



Surface tension (t_i can be obtained from the mem- 
brane area conservation law in the main approximation, 
A2(m) = 0. Formally, uq in (|A35p is order of 1, that is 
reflected by low index "0" . Note, that rotational part of 
external flow lu does not entered in (|A35p . This is the 
consequence of the fact, that the rotational part of the 
external flow enters in the exact equation on u though 
extension of time derivative, dt —* {dt — ujd^)^ that is 
corrections from oj to (jA35p arc formally order of a/A- 



4. Corrections order of s\^. 

In this subsection we find the first corrections in \/A 
to the equation of motion (jA35p . One should repeat the 
steps (|A24IA26IA28[ IA29IA32IA33p . accounting next or- 
der in -v/A and find correction iii order of -v/A to ii (jA35p . 

In the main approximation we required (|A24p . that 

[^17*]°"* = [^*u*]q" I ^ on the sphere with radius To . Here 
low index "0" stands for contribution from Xq,..., see 
definition (K3\i . whereas low index "1" in (|M8)) corre- 
sponds to contribution from Xi, . . .. Quantity Pv"^ calcu- 
lated on the membrane surface through found coefficients 
Xq, . . ., differs from that on the sphere: 



n=0 



5- [VvX 



(A37) 



where summation index n corresponds to contribution 
order of A"/^. Thus continuity of normal velocity com- 
ponent Vv^ in first order in \fK looks like 



6^ 



(A38) 



where indices "s" and "m" stands for sphere surface and 
membrane surface correspondingly. Condition can be 
rewritten as 



since = Xi and 

L J I m L J 1 1 



(A39) 



Next boundary condition Pl'^dkV^ = in the first order 
of \/A reads 



0, 



Yl'' + [Pl'^dkv'y^ = (A40) 

111 

where 
and 

Analog of (|A29p now is 

ui^-oVo^W^l+X^. (A41) 
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In (|A4ip velocity 



Full dynamical equation acquires the form 



(A42) 



where mq is from (|A35l) . 

Corrections to boundary conditions (|A32IA33[) are 

a-iH2 + (cto + o-o)-ffi + (TiHq = 

= LPiJL+L^iPoJL,, (A43) 

Vtai = [77 [P§,„] J 1^ + [77 [P§„] J 1^. (A44) 
In (|A43)l 

In (|I34)) 

h[/'§«a]J = (A45) 

S 

and 



lOuViJi 



+ (A-1) Au'P [0-i?i„] - 2u>o 



(A46) 



We do not write out boundary conditions on ^['^•°"' ^ since 
Z can be excluded from (jA44IA45P by taking divergence. 
On the way one obtains, that alternating part of surface 
tension 



^1- 0-1 1 [77 [r§,„]J'" 



r 



Now from boundary condition (jA32|) one can extract cor- 
rection ill to ii order of \fK: 



{dt-u:d^)u=^V^u + C{u) 



lO/i 
"a(2) 



+ MhU. 



(A48) 



5. Restriction onto sector / = 2. 

It follows from evolution equation, written in the main 
approximation over v^, that at long times all excess area 
becomes confined in harmonics 1 = 2. Coupling of higher 
order harmonics with sector / = 2 occurs only due to 
nonlinear terms, concerned in Subsection lA 41 Relative 
part of excess area, confined in higher order harmonics 
is order of A. Hence, back influence of higher order har- 
monics on sector I = 2 appears, if one accounts terms 
order of sA, that is the influence is negligible in our ap- 
proximation. Thus, to obtain corrections to dynamics 
in sector Z = 2 in the approximation, it is sufficient to 
account only coupling of mode / = 2 with itself. 

Bilinear function C{u) becomes in the case 



24(152 + 63A) ^ 
"^^"^ ^ (32 + 23A)^ " 



(A49) 



and nonzero elements of linear operator in the basis ((23 
are 



Ml 



1 _ 80^/3(56 - 31A) _^i5 _ 400(56 -41A) 



7(32 + 23A)2 



7V3(32 + 23A)2 ' 



7(32 + 23A)2 



ill = i^d^ + aa{0 + 2)u + a^iC{u) + MhU. (A47) 
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